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Abstract 

Within the framework of the recently proposed formalism using non-hermitean Hamil- 
tonians constrained merely by their VT invariance we describe a new exactly solvable 
family of the harmonic-oscillator-like potentials with non-equidistant spectrum. 
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1 Introduction 



An increase of interest in the VT symmetric Hamiltonians [Q - || may be explained 
by a number of their appealing properties. For illustration, let us pick up the most 
ordinary harmonic-oscillator Schrddinger equation in one dimension 



and change its variable r to x = r + ic, c > in such a way that x becomes treated as 
real, x G (—00, 00). Obviously, as long as r 2 = x 2 — 2icx — c 2 , the asymptotic growth 
or decrease of the old general solution ip(r) remains equivalent to the asymptotic 
growth or decrease of the new, complex function tp(x) = ip(x — ic). In the new, 
shifted bound-state problem 



boundary conditions remain unchanged, therefore. As a consequence, the spectrum 
E = E m = 2m + 1, m — 0, 1, ... of energies remains discrete, real and bounded from 
below, shifted merely by a constant c 2 > in the latter case. 

Our illustrative non-hermitean Hamiltonian commutes with the product of parity 
V (changing x to —x) and time reversal T (changing, formally, the imaginary unit 
i to —i). Daniel Bessis and Carl Bender with Stefan Boettcher conjectured 
that such a VT symmetry might be responsible for the reality of spectra for a much 
broader class of Hamiltonians. Their conjecture is widely supported by a number 
of tests. For the anharmonic V(x) = x 2 + igx 3 , its perturbative confirmation (viz., 
the proof using the Borel summability) has even been made available in the early 
eighties [§]. Very recently, the further numerical and quasi-classical evidence has 
been provided by the non-polynomial V(x) = x 2 {ix) 5 and by its supersymmetric 
partners 0J] as well as by certain hyperbolic [H and trigonometric M models. An 
additional, purely non-numerical backing of the hypothesis may, last but not least, 
rely upon the exactly solvable VT symmetric V(x) = exp(ix) jp and upon the quasi- 
exactly solvable polynomial V(x) = — x 4 + tax 3 + bx 2 + icx [0 and non-polynomial 
V(x) = x 2 + iax + (b + icx) /(l + dx 2 ) with d > [BJ . In the light of these and further 



- — + r 2 ) 4>(r) = Eip(r), ip(r) e L 2 (-oo, 00) 



dx 2 




<p(x) e L 2 (— 00, 00) 
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references |9|] it is rather surprising that no attention has been paid, up to now, 
to harmonic oscillators with a properly regularized centrifugal-like core of strength 
G = a 2 - 1/4 with a > 0, 

( d 2 G \ 

+ x 2 - 2icx + —- (p(x) = (E+c 2 ) <p(x), ip(x) e L 2 (—oo, oo). (1) 



dx 2 (x — ic 



The gap is to be filled by the present note. We are persuaded that this "forgotten" 
VT symmetric model deserves a few explicit comments at least. 



2 Quasi-parity 

In the usual one-dimensional and hermitean world the quadratic singularity proves 
too strong and, at nonzero G > 0, it cuts the real axis in the two separate, non- 



communicating halves ||10|| . In the present more permissive context, the threat is 
easily avoided by a shift of the singularity off the integration path. In the complex 
plane of x cut, say, from x = ic to x — > -Moo, our problem (U) remains well defined. 

What is equally important is its exact solvability in terms of the confluent hy- 
pergeometric special functions, 

ip(x) = C + (x- ic)~ a +i/2 e -(*-«) a /2 lFl (( 2 -2a- £)/4, 1 - a; (x - ic) 2 ) + 

+C_ (x - a C )"+i/2 e -(*-ic)72 iFi ^ 2 + 2a _ E y^ i + a; ( x _ ic f^ . 

For large \x\ this expression grows as exp(x 2 /2) and violates boundary conditions 
unless it degenerates to a polynomial (cf., mutatis mutandis, (TTJ). In this way we 
get the (complete) spectrum of energies 

E = E qn = An + 2 + 2qa 

numbered by the quasi-parity q — ± and integers n = 0,1,2,.... The related 
normalizable wave functions 

(p(x) = Const. (X - zc )-^+l/2 e -(^-i C ) 2 /2 L (-qa) _ ^ ^ 
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are defined in terms of the well known orthogonal Laguerre polynomials, 

LS(z) = 1, 
L{{z) =P + 1-z, 
L%(z) = (f3 + 2- z y-({3 + 2), 
L§(z) = (p + 3-z) 3 -3(f3 + 3)([3 + 3-z) + 2{/3 + 3), 

In the limit a — > 1 /2 and c — > our Hamiltonian re-acquires its hermiticity. Our set 
of solutions coincides with the well known one-dimensional harmonic oscillators and 
the quasi-parity degenerates to the ordinary parity, 

Vijj(r) = ip{—r) = qip(r). 

The spectrum of energies becomes equidistant, E +0 = 1, E_ = 3, E+i = 5, E-\ = 7 
etc. Precisely 2n+ (1 — ql)/2 real nodal zeros appear in the corresponding real wave 
functions. 

After we switch on a "subcritical", permitted central attraction on, the nodal 
zeros of f(x) in eq. (Q) will move upwards in the complex plane. Within interval 
— 1/4<G<0 with < a < 1/2, all the even and odd energies undergo an upward 
and downward constant shift, respectively. At the infinitesimally small extreme 
values of a ~ all the energies almost degenerate in doublets E± ~ 2, E±\ ~ 6, 
E± 2 w 10 etc. 

3 Strong repulsion and level crossing 

We have seen that the VT symmetry acts (or at least might act) as a simple and 
efficient means of regularization of a singularity in V(x) in one dimension ||. Our 
present example with a > 1/2 extends this idea in a way inspired by the solvability 
of the radial equation in three dimensions. There, G = l{t + 1) contains the angular 
momentum i = 0, 1, . . . and may become quite large. After the present regularization 
and analytic continuation of this model to one dimension the only important novelty 
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is the sudden disappearance of the (now, redundant) boundary condition in the 
origin. 

In the language using the complex coordinates x G (C we may also return from 
one dimension with x G (—00,00) to three or more dimensions with x G (0, 00). 
Beyond such a purely kinetic interpretation of our strongly repulsive core a genuine 
dynamical meaning of G = £(£ + 1) ^> 1 may be encountered, say, in nuclear physics 
where one has to use I = 3, 33, 117, 352, 517 and 1083 in an efficient approximative 
description of the respective nuclei A He, 16 0, 40 Ca, 90 Zr, 120 Sn and 208 Pb in the so 



called breathing mode [12 



After we turn on an enhanced repulsion in eq. ([!]) we discover a quasi- degeneracy 
and crossing of levels with opposite quasi-parities in the vicinity of every integer 
a = 1, 2, . . .. In the very first case with a = 1 we may factor the square (x — ic) 2 out 
of the states with the even quasi-parity, 

[(* - *c) 2 ] = 1, 
4~ 1} \(x-ic) 2 ] = -(x-ic) 2 , 
L { f 1] [(x - ic) 2 ] = -2(x - ic) 2 + (x - ic)\ 
L ( f 1] [{x - ic) 2 ] = -{x - ic) 2 [{x - icf - 6(x - ic) 2 + 6] 



The resulting formula [(x — ic) 2 ] = —(x — ic) 2 [(x — ic) 2 ] implies that the 
even and odd quasi-parity partners will coincide precisely at the "exceptional" JK| 
value of G = 3/4. An unavoided crossing of the energy levels occurs without their 
degeneracy. Similar phenomenon may be observed at all the subsequent integers 
a = 2, 3, . . . . 

During the steady growth of the repulsion a > 1/2 the relative displacement of 
the two halves of the spectrum distinguished by their quasi-parity is accompanied by 
certain interesting changes in the structure and position of the nodal zeros in <p(x). 
Their detailed analysis already lies out of the scope of the present note. Mathe- 
matically, it reflects a complex generalization of the usual Sturm-Liouville oscillation 
theorems Jl3| . 
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4 Weak core as a perturbation 

Besides the natural interpretation of small deviations from equidistant spectrum in 
a weak-coupling regime with G ~ we may also try to trim or suppress the influence 
of the core in eq. (|1|) via a sufficiently large screening c ^> 1. In such an alterna- 
tive setting our potential may be decomposed into its dominant (shifted) harmonic 
oscillator part V^ HO \x) — (x — ic) 2 and a well-behaved O(l/o 2 ) perturbation, 

V(x) = V {HO \x) + GW(x), W(x) = W {I \x) + W {II \x) + W {III \x). 

After a re-parameterization ji — g — c~ 2 and A = — c~ 4 the first, asymptotically 
dominant O (x~ 2 ) component of the anharmonicity 



x 2 + c 2 1 + g x 2 



appears quasi-exactly solvable at certain strengths G [LJ] . This correction has already 



been used in numerous methodical considerations [ITH]. The subsequent term 



is less common. It does not commute with the parity V and breaks the hermiticity 
of the (unshifted) oscillator, obeying only the overall VT invariance. The last, real 
and even component 

^ (m, M = -(^ = (V* 4 ) 

converts the perturbation W(x) to its present exactly solvable form. It is bounded 
and asymptotically decreasing. Its routine treatment, say, within the Rayleigh- 



Schrodinger perturbation formalism may be expected nicely convergent [16 . 



5 



5 Summary 



Many potentials of phenomenological interest are analytic functions. This makes 
(or at least might make) the underlying differential Schrodinger equation and many 
properties of its solutions much more transparent. In particular, we may imagine that 
all the functions which satisfy the equation on a real interval may be immediately 
continued into a bigger complex domain. 

In principle, the related possible shift or deformation of the axis of coordinates 
breaks the hermiticity of the Hamiltonian. Bound states acquire Im E ^ and be- 



come re-interpreted as unstable resonances [17]]. With quite a few important excep- 
tions: It is already known for many years that certain non-hermitean Hamiltonians 
H 7^ H + still do support perfectly stable bound states. The puzzling existence of 
these exceptional "stable resonances" with Im E = could prove helpful in phe- 
nomenological considerations and has been subject to an intensive study recently. 

In this context, our present note has shown that in the particular quantization 
scheme which weakens the hermiticity of a Hamiltonian to its mere VT invariance 
the one- dimensional superposition V(r) = r 2 + G/r 2 of the harmonic and centrifugal- 
like forces may be regularized by a purely imaginary shift of r in such a way that 
the whole model remains exactly solvable. 
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